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We demonstrate that the largest measured mass of a neutron star establishes an upper bound to 
the energy density of observable cold matter. An equation of state- independent expression satisfied 
by both normal neutron stars and self-bound quark matter stars is derived for the largest energy 
density inside stars as a function their masses. The largest observed mass sets the lowest upper limit 
to the density. Implications from existing and future neutron star mass measurements are discussed. 
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The number of neutron stars with measured masses 
has grown in recent years Q, ■ The most accurately 
measured masses are from timing observations of radio 
binary pulsars Q and, until recently, were consistent with 
neutron star masses in the range 1.26 to 1.45 M Q pj. 
Recent data on binaries containing pulsars and white 
dwarfs, however, indicate a larger range of masses 0- 
An example is the binary containing PSR J0751-1807, 
with 2.2 ± 0.2 M Q with la errors |4j. Data from x-ray 
binaries Q also suggest a wide range in masses, but are 
subject to greater theoretical and observational uncer- 
tainties. As neutron stars are expected to contain the 
densest cold matter outside black holes, the maximum 
neutron mass mass and the corresponding maximum en- 
ergy density are of great interest. 

We demonstrate here that a precisely measured neu- 
tron star mass sets an upper limit to the mass density, 
or, equivalently, the energy density, inside the star. The 
larger the measured mass, the smaller the density limit. 
A sufficiently large mass could delimit classes of possible 
equations of state (EOS's). A limit for this maximum 
density is preferred utilizing an analytic solution of Ein- 
stein's equations. This limit is checked by comparing 
numerical results for a variety of EOS's of both normal 
and self-bound stars. Recent neutron star mass measure- 
ments are summarized and inferences drawn. 

From the general relativistic structure equations 0, , 
the maximum compactness of a star is set by the limit 
R > (9/4)GM/c 2 0, where R and M are the stellar 
radius and mass, respectively. With the additional re- 
quirements that (i) nowhere in the star is the speed of 
sound c s greater than the speed of light c, (ii) c s is ev- 
erywhere real, and (iii) the EOS matches smoothly to 
calculable low density EOS's near the nuclear saturation 
density p s ~ 2.6 x 10 14 g cm" 3 , Ref. [8j showed that the 
compactness limit is increased to 



R > 2.MGM/c 



(1) 



This result improves the limit R <; 3.05GM/c 2 estab- 
lished 0| using t he p rescription c s = c above a fiducial 
energy density p / [lOj ■ The maximum mass inferred from 

— 1/2 

this prescription is proportional to p^ , but the com- 



pactness limit is independent of pf for pj << p Cl where 
p c is the central density of the star 0. 

The central mass density of a star must be greater 
than the average density p* — 3-M/(47ri? 3 ), the value for 
a uniform density star with the same mass and radius. 
Combining the compactness limit, Eq. l|T]l. with the con- 
stant density relation — p Cl yields 



p c .* ~ 5.80 x 10 15 (Mq/M) 2 



g cm 



(2) 



This is a plausible approximate lower limit to the central 
density p c for a star of a given mass, but it is not an 
absolute lower limit. (This lower limit cannot be made 
firm as causality has been imposed on a uniform density 
fluid in which transmission of signals is instantaneous.) A 
firm lower limit can be established, however, if an upper 
limit to R exists. One observational limitation originates 
from the most rapidly spinning pulsar, PSR B1937+21 
[Tlj . which has a fre quen cy v = 641 Hz. This leads to a 
lower limit to Mj R? [l2j and a lower limit 



p c , rot ~ 1.79 x 10 14 O/641 Hz) 2 g cm 



-3 



(3) 



which is, however, not very restrictive. A far more strin- 
gent limit could be achieved from a redshift observed 
from a neutron star. The largest observed redshift z b s 
sets a lower limit to M/R, implying 



Pc,z > 



^obs 

(2 + z obs ) 
AttM 2 V 2G(1 + Zobs y 



(4) 



Recently, z \, s — 0.35 was reported [13J for the x-ray 
bursting source XTE J1814-338. With this value, 



p c , z > 1.69 x 10 15 (M /M) 2 g cm 



-3 



(5) 



The central question is, how much greater can p c be 
compared to any of the above expressions for physically 
motivated EOS's? If this question can be answered, an 
upper limit to the density inside a star of a given mass 
can found. An important consequence of the existence 
of an upper limit is that the largest measured neutron 
star mass would set an upper limit to the density of cold 
matter. (In a dynamical environment, such as the gravi- 
tational collapse of a stellar core to a black hole or a high 



2 



energy heavy ion collision, matter becomes hot and may 
achieve higher densities.) An additional consequence is 
that one could infer whether or not non-nucleonic de- 
grees of freedom, such as hyperons, Bose condensates or 
quarks, which generally reduce the maximum mass, can 
exist in the cores of neutron stars. 

Some insights can be gained by comparing analytical 
solutions to Einstein's equations with numerical solutions 
employing model EOS's. The known analytic solutions 
fall into two classes: (i) the class that describes "normal" 
neutron stars for which p c vanishes at the surface where 
the pressure vanishes, and (ii) the class that describes 
"self-bound" stars for which p c is finite at the surface. 
In the first class, there are only three known analytic 
solutions: the Tolman VII solution la, B uchdahl's solu- 
tion 01, and the Nariai IV solution |l5|. In the second 
class, an infinite number of analytic solutions exist, but 
the useful ones are variants of the Tolman IV and VII 
solutions 0, 0] , as well as the uniform density case 0] . 

All known analytic solutions are scale-free; they 
depend parametrically on the compactness ratio 
(3 = GM/Rc 2 . However, by coupling these solu- 
tions with Eq. i.e., by setting the compactness 
(3 = j3 c = 1/2.94, one can obtain relations between p c 
and M analogous to that of a uniform density fluid, Eq. 
(0. We thus need to relate the central density to f3 and 
M for these solutions. These are summarized below: 

1. Tolman VII: This solution stems from the ansatz [f| 



P = Pc [1 - (1 - w)(r/R) 2 } , 



(0) 



where the parameter w — p s j p c is the ratio of the energy 
densities at the surface and the center, which can vary 
between and 1. (The case w = 1 represents the uniform 
density fluid.) This leads to 

j .) / <-. i, 

Pc,VII 



4tt(2 + 3w) V G 
1.45 x 10 16 fM c . x 1 
(1 + 1.5iu) 



1 

AI 2 



L 

M 



g cm 



(7) 



This solution is valid, in the case w = (the normal 
neutron star case), for (3 < 0.3862 ~ 1/2.59. For positive 
w (the self-bound case), the solution is valid for larger 
values of (3. Thus, this solution is useful for the case 
(3 = (3 C . For w > 0, the central density decreases for a 
given mass relative to the normal neutron star case. 



2. Buchdahl: This solution uses the EOS ansatz pc 2 = 
\2\Jy*P — 5P, where is a constant. In this case 

n\;i ■ ■ :> VMJ - ),-)- / <" ), 1 " 

Pc.Buch 



1 - Wc 



G 



3.89 x 10 15 (Mq/MY g cm 



1 

M2 

(8) 



However, because Buchdahl's solution is invalid when 
/3 > 1/5, the value for which the central sound speed 



becomes infinite, it cannot be used for the case [3 = j3 c . 



3. Nariai IV: This solution is characterized by [l 

Pc,NIV 



3 

8^ 



(6 ex.') 

(a — 1) cos \/3f3 c H sin \/3f3 c 



G 



1 

M 2 



9.88 x 10 15 {Mq/MY g cm- 



Here, 



a 



12 + /? c 



2 + f3 c + 2y/l - 2fi c 



(9) 



(10) 



This solution is valid for (3 < 0.4126, the value for 
which the central pressure and sound speed become 
infinite. This solution can also be generalized to include 
self-bound stars, and as for the Tolman VII case, the 
central density for a given mass decreases from that of 
Eq. © as the ratio w is increased from 0. 

4- Tolman IV (generalized): Lake |l6j showed that the 
ansatz for the metric function 



,u{r) _ 



[N - (3(2N + 1) + (3(r/R) 
N N (l-2f3) N ~ 



2lJV 



(11) 



where N is a positive integer, produces an infinite family 
of analytic solutions of the self-bound type. Four of these 
were previously known (N = 1,3,4 and 5). The case 
N = 1 cannot properly be applied to our problem as this 
solution is finite only for (3 < 1/3. The most relevant case 
is for N = 2, for which c s sa y/l/Z throughout most of 
the star, similar to the behavior of strange quark matter. 
For this case, 



Pc.TIV — 



3 

4tt 



2-2/3. 
2~5f3 c 

116 



2/3 



G 



1 

M 2 



1.56 x 10 ib (Mq/MY g cm" 3 . 



(12) 



This solution is valid for (3 < 2/5. The ratio of the surface 
to the central densities for the case N = 2 is 



6-10/3 ^2-5/3 
2-2/3 



2/3 



(13) 



which is approximately 0.32 for (3 = (3 C . With increasing 
N, the central density for a given mass decreases from 
that of Eq. {T5J|, and in the limit N » 1, p c ,Tiv(M) 
approaches the uniform density result. 

To investigate the relevance of analytic relations be- 
tween the central density and mass, we carried out 
numerical integrations of the TOV structure equations 
for a multitude of EOS's, including potential and field- 
theoretical models, and models that contain strong soft- 
ening due to the occurence of hyperons, Bose condensates 
or quark matter, including the case of self-bound strange 
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FIG. 1: The central energy density and mass of maximum 
mass configurations. Symbols reflect the nature of the EOS's 
selected from Refs. Ml. NR are non-relativistic potential 
models, R are field- theoretical models, and Exotica refers to 
NR or R models in which strong softening occurs, due to the 
occurence of hyperons, a Bose condensate, or quark matter. 
The Exotica points include self-bound strange quark matter 
stars. For comparison, the central density - maximum mass 
relations for the redshift, Tolman VII [w = 0), and Tolman 
IV (N = 2) bounds from Eqs. ©, Q, and 1121 are shown. 
The dashed lines for 1.44 and 2.2 Mq serve to guide the eye. 
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quark matter stars. The EOS's were chosen from Refs. 
18]. Fig. ^displays the maximum masses and the central 
energy densities of the maximum mass configurations. 
Analytic solutions for the Tolman VII normal neutron 
star case (w = 0), the Tolman IV case for N = 2, and 
the uniform density case are also displayed. The paths 
of other analytic solutions are not shown for clarity; they 
scale up or down from the analytic results shown (all 
analytic solutions for fixed obey the law p c oc M~ 2 ). 

It is fortuitous but significant that the Tolman VII 
solution forms a strict upper limit to the density of a 
maximum mass star, for each of the EOS's displayed. 
We therefore conjecture that the Tolman VII curve marks 
the upper limit to the energy density inside a star of the 
indicated mass. Since the maximum density achieveable 
with a given EOS is the central density of the maximum 
mass star, a stellar mass measurement can be directly 
converted into an upper limit for the maximum density. 
Since the measured mass must necessarily be less than 
the neutron star maximum mass, this limit also forms 
an absolute upper density inside any compact star. In 
other words, except in transient situations such as the 
Big Bang or in relativistic heavy ion collisions, this curve 
displays the ultimate density of matter. 

The above results are given in terms of the central mass 
or energy density p. However, most models of dense mat- 
ter are formulated in terms of the baryon number density 
n. Therefore, it is of interest to examine relationships be- 



FIG. 2: Measured and estimated masses of neutron stars in 
radio binary pulsars and in x-ray accreting binaries. Sources 
are listed by letter (Refs. 0, Q and [20-35]). Error bars 
are la. Vertical dotted lines show average masses of each 
group (1.62 M©, 1.34 M© and 1.56 Mq); dashed vertical lines 
indicate inverse error weighted average masses (1.48 Mq, 1.41 
M and 1.34 M Q ). 

tween the central baryon density and mass for the max- 
imum mass configurations. A good rule of thumb for 
converting p to n, using n s ~ 0.16 fm~ 3 , is 

p/ps « 0.9(n/n a )[l + 0.11(n/n s ) 3 / 4 ] . (14) 

The number density so obtained is indicated on the top 
scale of Fig.^ We emphasize that the plotted points are 
positioned using p c , not n c , in this figure. 

The most accurately measured neutron star masses are 
from timing observations of radio binary pulsars 0, • 
These include pulsars orbiting another neutron star, a 
white dwarf or a main-sequence star. 

Measured masses are summarized in Fig. [21 and are 
plotted with la uncertainties. Ordinarily, observations 
of pulsars in binaries yield orbital sizes and periods from 
Dopplcr phenomenon, from which the total mass of the 
binary can be deduced. But the compact nature of sev- 
eral binary pulsars permits detection of relativistic ef- 
fects, such as Shapiro delay or orbit shrinkage due to 
gravitational radiation reaction, which constrains the in- 
clination angle and permits measurement of each mass 
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in the binary. The largest accurately measured mass 
originates from the binary pulsar system PSR 1913+16, 
whose masses are 1.3867 ± 0.0002 and 1.4414 ± 0.0002 
M Q , respectively [34^. 

A significant development concerns mass determina- 
tions in binaries with white dwarf companions, which 
show a broader mass range than binary pulsars having 
neutron star companions. Ref. suggests that a nar- 
row set of evolutionary circumstances conspire to form 
double neutron star binaries, leading to a restricted range 
of neutron star masses. The implication of this restriction 
for other binaries remains to be explored. A few cases of 
white dwarf binaries that contain neutron stars consid- 
erably larger than the canonical 1.4 M© value have been 
reported. A striking case is PSR J0751+1807 in which 
the estimated mass with la error bars is 2.2 ± 0.2 M©. 
For this neutron star, a mass of 1.4 M© is 4er away. If this 
mass determination holds up after further observations, 
the central density constraints become intriguingly close 
to the estimated density for the quark-hadron phase tran- 
sition. Raising the limit for the neutron star maximum 
mass could also mark the boundaries of other families of 
EOS's in which substantial softening begins around 2 to 
3n s . This is significant, since exotica generally reduce 
the maximum mass appreciably. 

The simple mean of the measured neutron star masses 
in white dwarf-neutron star binaries exceeds that of the 
double neutron star binaries by about 0.22 M Q (Fig. |5J). 
Nevertheless, caution is in order since the 2a errors of 
most of these systems extend below 1.4 M©. Continued 
observations are required to reduce these errors. 

Masses can also be estimated for binaries which con- 
tain an accreting neutron star emitting x-rays. Some 
of these stars are characterized by relatively large masses 
but also large estimated errors (Fig. EJ). The system Vela 
X-l is noteworthy, because its lower mass limit (1.6 to 1.7 
M©) is constrained, albeit mildly, by geometry j^. The 
source 4U 1700-37 might be a black hole, due to lack of 
oscillations in its x-ray spectrum |20|. Another object, 2S 
0921-630 could either be a high-mass neutron star or 
a low-mass black hole. Although not yet demonstrated, 
it is widely believed that black holes formed in gravita- 
tional collapse have masses that exceed the neutron star 
maximum mass. These latter two objects could play a 
significant role in determining the neutron star maximum 
and the black hole minimum masses. 
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